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We present a detailed quantum mechanical treatment of the photodetachment of 
H by a short laser pulse in the presence of crossed static electric and magnetic fields. 
An exact analytic formula is presented for the final state electron wave function (de- 
scribing an electron in both static electric and magnetic fields and a short laser 
pulse of arbitrary intensity). In the limit of a weak laser pulse, final state electron 
wave packet motion is examined and related to the closed classical electron orbits in 
crossed static fields predicted by Peters and Delos [Phys. Rev. A 47, 3020 (1993)]. 
Owing to these closed orbit trajectories, we show that the detachment probability 
can be modulated, depending on the time delay between two laser pulses and their 
relative phase, thereby providing a means to partially control the photodetachment 
process. In the limit of a long, weak pulse (i.e., a monochromatic radiation field) our 
results reduce to those of others; however, for this case we analyze the photodetach- 
ment cross section numerically over a much larger range of electron kinetic energy 
(i.e., up to 500 cm -1 ) than in previous studies and relate the detailed structures 
both analytically and numerically to the above-mentioned, closed classical periodic 
orbits. 
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I. INTRODUCTION 



High resolution studies of atomic Rydberg states in the presence of external static electric 
and magnetic fields have proved to be exceedingly fruitful for the investigation of atomic 
dynamics because, owing to the large radial extent and weak binding of atomic Rydberg 
levels, the effects of external static fields are much more significant for Rydberg levels than for 
atomic ground or low- lying excited states Q]. Consequently for more than a quarter century 
(up to the present) experimentalists and theorists have been investigating atomic Rydberg 
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signatures in Rydberg spectra |2l|, |2j|, of doubly-excited states in crossed fields [22! , °f 
recurrence spectra [sO, Lul [3] and of various aspects of electron dynamics in combined 



Coulomb and crossed static electric and magnetic fields [Jj, LUj [la [la, [121, 12 
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The related problem of photodetachment of a weakly bound electron (e.g., as in pho- 
todetachment of a negative ion) in the presence of crossed static electric and magnetic fields 
has been the subject of fewer investigations despite its having a comparably rich spectrum. 
(Note that the weakly bound electron in a negative ion can simply decay, or become de- 
tached, solely due to the presence of the external static electric and magnetic fields, a process 
that has long been studied theoretically, as in, e.g., J37], 38].) Experimentally, crossed field 
effects have been found to be significant in photodetachment of negative ions in the presence 
of a static magnetic field owing to the influence of the motional electric field experienced 
by the detached electron HQ. The photodetachtnent spectrum of H" in the „ of 
crossed static electric and magnetic fields has been treated theoretically by Fabrikant [41J 
and by Peters and Delos 3, Q) ; a generalization to the case of photodetachment of H~ in 



the presence of static electric and magnetic fields of arbitrary orientation has been given by 
Liu et al. 44,Q,Q- 

In each of these works the static fields are assumed to be sufficiently 
weak that they do not affect the relatively compact initial state. Fabrikant gave the 
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first quantum treatment of single photon detachment in crossed static electric and magnetic 
fields using the zero-range potential model to describe the initial state of H~; rescattering 
of the electron from the potential was also investigated, although the effect was found to be 
small except for high magnetic field strengths. Peters and Delos [3] gave a semiclassical 
analysis of H~ photo detachment in crossed fields and correlated significant features of the 
spectrum with closed classical orbits. Subsequently they derived quantum formulas for this 
process (using the zero-range potential model for the initial state) and exhibited the con- 
nection to their predicted classical closed periodic orbits j^j]. The generalization of Liu et 



al. 



4J,|45 



46l | to the case of static electric and magnetic fields of arbitrary orientation is also 



based upon the zero-range potential model. In all of these works the electromagnetic field 
that causes photodetachment is assumed to be weak and monochromatic. Also, the pho- 
todetachment spectrum is analyzed numerically only over a very small energy range above 
threshold. 

In this paper we consider detachment of H~ by a short laser pulse in the presence of 
crossed static electric and magnetic fields. We present an analytic expression for the final 
state of the detached electron taking into account exactly the effects of the laser field as well 
as both static fields. The initial state is described by the solution of the zero-range potential, 



as in all other quantum treatments to date 



41 



43 



45|,|46!. We present also an analytic 



expression for the photodetachment transition amplitude that can be used to describe the 
probabilities of multiphoton detachment in crossed fields. In this paper, however, our focus 
is on single photon detachment by short laser pulses and on the connection between the 
detached electron wave packet motion and the predicted classical closed periodic orbits of 



Peters and Delos 



42|. As noted by Alber and Zoller 



47| (in connection with electronic 



wave packets in Rydberg atoms), such wave packets "provide a bridge between quantum 
mechanics and the classical concept of the trajectory of a particle" and "the evolution 
of these wave packets provides real-time observations of atomic or molecular dynamics." 
We show this connection for the case of short pulse laser- detached electron wave packets 
in crossed static electric and magnetic fields. In addition, we show analytically how our 



4l| 



short pulse results reduce to the quantum monochromatic field results of Fabrikant 
and Peters and Delos j^J in the long pulse limit as well as the connection between our 
analytic quantum formulation for the photodetachment spectrum and those features that 
we associate with the predicted classical closed orbits [42 1. Finally, we present numerical 



i 



results in the long pulse limit over a large energy range above the single photon detachment 
threshold in order to demonstrate clearly these manifestations of classical behavior in our 
predicted photodetachment spectrum. 

This paper is organized as follows: In Sec. II we present our theoretical formulation 
for detachment of H~ by a short laser pulse in the presence of crossed static electric and 
magnetic fields. In particular, in this section (with details given in an Appendix) we present 
an exact, analytic expression for the wave function for an electron interacting with both the 
laser pulse and the crossed static electric and magnetic fields. We present here also analytic 
expressions for the transition probability amplitudes for both a single laser pulse and a double 
laser pulse (i.e., two coherent single pulses separated by a time delay). In addition, the long 
pulse (monochromatic field) limit of our results is presented and this result is compared 
with a number of prior works for various static field cases. In Sec. Ill we establish the 
connection between the long pulse limit of our results and the closed classical periodic orbits 
predicted by Peters and Delos j3| . In Sec. IV we present our numerical results, starting first 
with a comparison with prior results for the long pulse (monochromatic field) case and then 
examining the short pulse case, including the final state motion of the detached electron 
wave packets. 



II. THEORETICAL FORMULATION 



We consider photodetachment of H~ by one or more short laser pulses in the presence 
of crossed static electric and magnetic fields. In the final state, we assume the detached 
electron experiences only the laser and static fields; we ignore final state interaction of the 
electron with the residual hydrogen atom. For weak external fields, this is expected to be 
a good approximation for this predominantly single photon process. In this section, we 
first give the S'-matrix transition amplitude for photodetachment of H~. Then we present 
an exact quantum mechanical solution to the time-dependent Schrodinger equation for the 
final state of the detached electron in both the crossed static electric and magnetic fields 
and the time dependent laser pulse. We then use this result together with S-matrix theory 
to obtain detachment rates and cross sections. Atomic units are used throughout this paper 
unless otherwise stated. 
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A. S'-matrix Transition Amplitude for Photodetachment of H 

We adopt the Keldysh approximation for the final-state, i.e., we neglect the binding 
potential |4£|. In this case, the S'-matrix transition amplitude from the initial state ipi to 
the final state tjjf is given by 

/oo 
-oo 

where Vj represents the laser-electron interaction and the bracket () stands for integration 
over momentum space. For the zero range potential for which the bound state wave function 
has the form e~ Kr /r, the S'-matrix element in Eq. can be shown to be gauge-invariant j3, 
0]. Such a bound state wave function can be used to represent the weakly bound electron 
of H~. We use that of Ohmura and Ohmura j^jj, which in momentum space is given by 



where Cj is a normalization constant and £j is the initial state energy. Using the variational 



results of Ref. [51| and effective range theory for a weakly bound s-electron |52j, one ob- 
tains (53! Ci = 0.31552 and = —0.027751 a.u.. The gauge-invariant S-matrix transition 
amplitude for H~ detachment is then given by (cf. Eq. (27) of Ref. j^|) 

(% W = 1 r dt ' (Mp^IQ/V^) e~^'. (3) 

J —CO 

B. The Final State Wave Function 

In order to calculate the S-matrix transition amplitude in Eq. Q, we present in this 
section an analytical expression for the final state wave function ipf. As aforementioned, 
we neglect the binding potential after detachment. Therefore, ipf is actually a Volkov-type 
wave function that describes a free electron moving in the combined field of the crossed 
static electric and magnetic fields and the time- dependent electric field associated with the 
short laser pulse. In Fig. ^we illustrate the configuration of the external fields in which the 
detached electron moves: the uniform static magnetic field defines the x axis and the static 
electric field defines the z axis, i.e., 

B = Bi (4) 
E s = E s k. (5) 
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FIG. 1: Geometrical arrangement of fields in photodetachment of H~ by a linearly polarized laser 
(with electric field El) in the presence of crossed static electric (Eg) and magnetic (B) fields. Both 
the laser and the static electric fields point along the z-axis. As indicated, the drift motion of the 
detached electron is along the y-axis. 

We assume that each laser pulse has the following general form, 

E L (t) = E e- a2(t - r)2 sin {ut + (3) k, (6) 

where u is the laser frequency, r is the time delay with respect to t = 0, and (3 is a 
(generally constant) phase. The duration of the laser pulse is defined to be the full width 
at half maximum (FWHM) of the laser intensity, and is given by 

T p = V2 h2/a. (7) 
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We introduce the vector potentials for the magnetic field and the laser field respectively, as 
follows: 



At 



■zBj, 



E h (t')dt' 



(8) 
(9) 



where c is the speed of light in vacuum. 

The final state wave function for the detached electron is obtained as the solution of the 
time-dependent Schrodinger equation (TDSE) in momentum space, 

9 +\ _ TJjh), 



i^f\p,t)=H^(p,t), 



(10) 



in which the Hamiltonian H is given by 

1 



H(p,t) 



p + -{A L + A B ) 

c 



r • E< 



1 2 9 2 . ( E s \ d 1 2 1 , , 



(12) 



where oj c = B/c is the cyclotron frequency. 

It can be shown that Eq. (Jl(Jj) has an exact analytical solution. The details of the 
derivation are presented in Appendix [X] The final expression of the solution is given by 



ipf\p,t) = 8{p x - k x )5{p y - ky)exp[-ie f t -if{t)] 



X 04, 



-1/4, 



V2L 



exp 



-ib(ky,t)V2( p 



in which g n _ (x) is defined by 

9n z (x) 



-x 2 /2rr 
±1 ri 



where H Uz (x) is the n z th Hermite polynomial. In Eq. ()13|) we have also defined 



(13) 



(14) 



2 

2^ 



/(*) = + ^ / A 2 L (t')dt'-[ L(t')dt', 

&(M) = C* - &)/vkf. 



(15) 

(16) 
(17) 
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where the arguments ( Pz and C,k y are given by 



Cp,(*) = Ip*-£(*)]/V2^ (18) 

( ky = 4= (ky ~ — J , (19) 



and the energy of the n^t/i Landau level e nz in Eq. (|15j) is given by 

n z + cu c . (20) 

In Eq. (116(1 . L(t) and £(t) are functions related to the vector potential of the short laser 
pulse. L(t) is defined as 

L(t) = ^\t) - l -e{t) - \A L {t)i{tl (21) 



while satisfies the following differential equation: 



2 



m+<4t(t) = --£-Mt), (22) 

where £(t) denotes the second derivative of £(£). We present the exact solution for in 
Appendix El However, in the long pulse case (a/u <C 1), a simplified expression for 
can be obtained as 

~ a (u) e- a ^ l - T)2 cos (cut + /?) , (23) 

where we have defined 



2 



«M= r^S - (24) 

cj (w 2 - ^ 2 ) 

In this case, the first derivative, is thus given by 

i(t) ~ -cua (w) e -° 2 (*- T ) 2 sin (cut + /5) . (25) 

In order to investigate wave packet dynamics, it is useful to derive an expression for the 
final state wave function in which its z component is given in coordinate space. This is 
achieved by taking the Fourier transform of Eq. (fT3*|) with respect to p z , i.e., 

i r°° 

4jf( Px , Py ,z,t) = 4% y ^(p x ,p y ,z,t) = — / dp^f(p,t)e 1 ^. (26) 

Changing the integration variable to ( Pz (cf. Eq. (|18(l). we obtain 

i>f(jPx,Pv,z,t) = S(p x - k x )5(p y - k y )i nz uj x J^g nz \^fuj~ c z - b(k yj t)} 

x exp [iz£(t) - is ft - i fit)} , (27) 

where we have made use of Eqs. 7.388(2) and 7.388(4) in Ref. 
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C. 5-matrix Amplitude for Photodetachment of H" 



In order to examine the motion of the detached electron wave packet in crossed E and B 
fields, we define in analogy to Eq. (23) of Ref. J3] a time-dependent transition amplitude 
Rfi(t) from the initial state to the final state (k x , k y ,n z ): 



V2vr 



: / uj r . 

-oo 



- x ^g nz (V2( Pz ) exp \ib{k y ,t')V2( p ] dC P 



(28) 



where eji — £/ — £j, b(k y ,t') is given by Eq. (fTTjl . and where we have used Ea. (TEfl) for the 
final state wave function in Eq. (j3J). Using Eqs. 7.388(2) and 7.388(4) in Ref. [5J| to carry 
out the integration over ( Pz , we obtain 



;n z +l 



a 



dt'e^'+^u^g^ [b(k y ,t')]. 



(29) 



Note that in the limit of t — > oo, Rfi(t) reduces to the 5"-matrix transition amplitude 



i.e.: 



(30) 



In principle, with this analytical S'-matrix amplitude one can readily calculate the total 
and multiphoton transition rates, as done for H~ detachment in a static electric field in 
Ref. 0|. However, in the present paper, we restrict our consideration to the one-photon 
detachment process (Note that there are still many cycles in the laser pulses that we will 
consider in this work). Our analytical results facilitate easy comparison with some other 
previous results. Consequently, we evaluate Eq. ()29J) only to first order in the laser electric 
field strength Eq, i.e., we employ the following approximations: 

e i/(f) ~ 1 +«'/(*') 



1+t— «(0. 

fU) r . 



9U Z [b(ky,lf] ~ guACky) + g'n z {(ky 



where g' n _(Ck y ) stands for the derivative of g nz (Ck y )- Thus, to first order in E , the time- 
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fix fiyTlz 



dependent transition amplitude is given by 

fix fcyflz 

J-oo V W c 



(31) 

Note that, as usual, the first term in Eq. (p?Tj) does not contribute to the photodetachment 
process (since for t — > oo, the only contributions are for Efi — > 0); hence this term is discarded 
in the following discussion. 

In the long pulse approximation, with the help of Eqs. (|2*3J) and (|25j). one can show that 



I ^iQnA^y) 5/4 



UJ, 



'i nz CiC y g nz (Cky 



! dt'e- a2{t '~ T)2+i£ ^' sin {ujt' + /3) 

J —oo 

l ^ f dt'e- a2 ^-^ +i£ ^cos(ujt' + P), (32) 

f J —oo 



UJc J —oo 



which reduces to 



Rfht) 



-i n *C, 



t—J/I M.(CO + u c ( y g nz {( ky )] 



2uj"r 



(33) 



if we neglect the emission process (i.e., if we discard terms involving e +luJt ' 
The integration over t' in Eq. (|3*3*|) can be carried out analytically: 

e ~a 2 (t' -t) 2 +ie fit 1 -iu>t' -ij3 



2a 
x exp 



erf 



a (t — t) — i 



2a 



uj 



4a 2 



+ i (e fi — uj) t — i(3 



Thus, for the single laser pulse in Eq. (JHJ), the first-order time-dependent transition amplitude 
in the long pulse approximation is given by 



2ul' A 



xD sg i(e fi ,t)5 a (e f i 



uj 



(34) 
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1 + erf 



where we have defined 

D sgl (e fl ,t) = e< £ f>-^ 
and have also introduced the quasi-5-function |5 

S a (sfi ~ w) — (27r 1 ^ 2 a) exp 



a (t 



.e fi - uj 

i— 

2a 



Aa 2 



(35) 



(36) 



In the limit that our finite laser pulse becomes a monochromatic plane wave, the quasi-5- 
function becomes the usual Dirac ^-function, 



5(e fi -u) = \im.5 a (e fi - uj 



(37) 



Taking the limit t — > +oo, we obtain from Eq. (J34"jl the following analytical expression 
for the S'-matrix amplitude for the case of a single, finite laser pulse: 

s si na(uj) 



c(l) 



LO, 



5/4 



xeV' i-w 



6 }i - UJ) 



(38) 



where we have used the fact that erf(oo + iy) = 1 for any finite real number y. 



D. Detached Electron Wave Packet 



We may obtain the detached electron wave packet probability amplitude as a sum over 
all final states of the product of the time-dependent transition amplitude for transition to 
the final state (k x , k y ,n z ) at a particular time t [Eq. ()34|)] and the wave function [Eq. (|27j)] 
for that state (cf. Sec. II. D of Ref. 

. ■ (39) 



3): 

°° roo poo 

lh/p(Px,Py,Z,t) = / dk * / d K^klk y ,n\Px^Py^,t)\R { l}{t) 

By using Eqs. (|27j) and ()34|). the wave packet for the single laser pulse (JSJ) is given by 



ipWp(Px,P y ,z,t) = -C ~ - ^ (-l) n 'exp 



2uJ r 



n z =0 



iz£ - is' f t - %——(;-£ 



X9n z - b{p y , t)} [ujg' nz (C Py ) + v c ( y g nz (( Py )] 

xD sgl (e' fl ,t)5 a (e' fl -uj) 



(40) 



where ( Py and b(p y , t) are defined by Eqs. (fT^j) and (fT7|) respectively, e'u = e'j- — and a'* is 
given by Eq. with (k y replaced by £ Py . 
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E. S-matrix and Wave Packet Amplitudes for the Double Pulse Case 



We consider here the case that there are two laser pulses of the form of Eq. pjl. with the 
second one delayed with respect to the first by a time interval r and having a relative phase 
of P, i.e., 



E£ bl (t) =E e~ a * sin (cut) + e 



-a 2 (t-r) 2 



sin {ojt + (3) 



k. 



(41) 



To first order in Eq, it is easy to show that for the double laser pulse case, the time-dependent 
transition amplitude is given by 



dbi 



4^a^[ug' nz (C ky )+u; c C y g nz (C ky )] 



1 2u^ 



x 



Ddbi( £ fi, t)5 a [Sfi 



to 



(42) 



where the function Z?dbi(£fi> t) is given by 



D d bi(£fi,t) = 1 + erf 



at — i 



2a 



! ft -J)T-ifl J 1 _|_ er f 



a (t — r) — i 



2a 



When £ — > oo, the above formula reduces to 

dbl / , 1X \ dbl 

I 

kx hyTiz 



/ k x k y n z 



lim 



8 a (e fi - u) . 



(43) 



(44) 



The wave packet amplitude for the double laser pulse case is correspondingly given by 



^(p x , Py ,z,t) = -C™^ £(-!)»- exp 



2u r 



n,=0 



izi - ie'ft - i^=( Py £ 



*9n z [V^~cZ ~ b(p y , t)] [wg' nt ((py) + WcCv9n M (£py)] 

xD m (e' fi ,t)S a (e' fi - u) . 



(45) 



F. Photodetachment Cross Section 



The transition probability to a particular final state (k x , k y ,n z ) is given by 



(46) 
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and the total photodetachment probability is calculated by integrating over all final states, 

00 roc poo 

P=Y. dk * dk yPk x k ynz . (47) 



For an infinitely long, monochromatic beam, the probability P is proportional to time, t. In 
this case, it does not make sense to talk about the total transition probability. Instead, one 
normally considers the total transition rate, W, which is given by 3] 

W = lim -P. (48) 

t— too t 

The total photodetachment cross section is obtained by dividing the total photodetachment 
rate W by the photon flux F (the number of photons per unit area per unit time): 

W 

&pw p 1 (49) 

where 'pit;' stands for the monochromatic plane wave case. 

For the short laser pulse case, it does not make sense to talk about a transition rate since 
the transition probability is not simply proportional to time t. In addition, the photon flux 
F is not well defined. Nevertheless, it is possible to renormalize the total probability for 
detachment by a short laser pulse in such a way that the renormalized probability reduces, 
in the limit of an infinitely long pulse, to the usual formula for the photodetachment cross 
section. Since the renormalized probability will have the dimensions of area, we denote 
it as an effective photodetachment cross section, a . To derive this effective cross section, 
one uses the time duration of the laser pulse as the unit of time. One calculates the total 
photodetachment probability P during the laser pulse duration and the total number of 
photons per unit area (i.e., the photon density), E, during the laser pulse duration. Then 
an effective photodetachment cross section, a, may be defined as 

a = |. (50) 

Clearly a defined in this way has the dimensions of a cross section. In the rest of this paper, 
a should be understood to be this effective photodetachment cross section, i.e., calculated 
according to Eq. (|50|). We shall show below that this a for the short laser pulse case 
reduces in the limit a — > (cf. Eq. (jBJ)) to the usual photodetachment cross section for a 
monochromatic plane wave. 
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It has been shown by Wang and Starace |55( that for a single Gaussian pulse defined by 
Eq. © with r = j3 = 0, the photon density £ is given by the following formula: 



cEl v 7 ^ 



8ttuj 2a 

For the double pulse case (cf. Eq. (jSj) ), £ is correspondingly given by 



^dbl 



CE 2 V 7 ^ ^ + COS/gexp (_ a 2 r 2 /2) ] 



(51) 



(52) 



57TOJ a 

Taking (3 and r to be zero in Eq. (|38|) . and using Eqs. (J24j) and (|46|) - (J51j) . we have for the 
photodetachment cross section of H~ by a single pulse of the form of Eq. ©: 

4tt 2 c:v 



a 



(i) 



Ojj(u! 2 — Ul 2 ) 2 



X 



n z =0 

dk x 5 a (e fi -Lo), 
where we have employed a second quasi-5 function 



/ d ^y [^c9n z {(k y )+ug' nz {Cky)]' 

„ —n J —oo 



S a (e fi - u) 



1 



exp 



2a 2 



(53) 



(54) 



«V 27T 

which reduces to the usual Dirac ^-function in the limit of a monochromatic plane wave, 
i.e., 

5 (Efi — to) = lim 5 a (Efi — uj) . (55) 

a — >0 

Note that 

e fi -u= l -[kl + Q(C kv )], (56) 
in which we have defined (cf. Eq. ITo^) 

Q ft*,) = ^ {(k y + Cmin) , (57) 

where 



(n z + -)uj c H 

v « t 2 ; 2c^ 



£i - 



(58) 



The integration over fc^ in Eq. (J53|) has an analytical result when Q (Ck y ) > 0. The result is 



dk'r 



V2a% 



exp 



: fx 



UJ 



V2o% j 



dk x exp 



2a 2 

(kl + Q(C ky )Y 



(59) 



1 



^5 



(Cfc„ + Cmin) exp 



(Cfc H + Cmin) 

4o; r a 2 



Ki 



Es {Cky + Cmin) 

4ccVQ; 2 
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where we have used the following formula (cf. Eq. (3.323) on p. 307 of Ref. |54j): 



tteexp [-P 2 x a - 2 7 V] = T-^-e^l^Kr ( 7 4 /2/3 2 ) 

P 



which holds for |arg/3 
(cf. p.375 of Ref. 
numerically. 



< j and |arg7| < ~, and where K u (z) is a modified Bessel function 
When Q (Ck y ) < 0, the integration in Eq. (J5T?|) must be done 



1. Plane Wave Limit of the Cross Section 



In the plane wave limit, a — > 0, the integration over k x (making use of Eq. (jM|)) becomes 



dkj {e fi -u) = J dkj (h 2 x + l -Q (C fc J J 



(60) 



This integral is non-zero only when Efi — uj = \k 2 x + \Q ((k y ) = 0, i.e., when we have strict 
energy conservation. For non-zero real k x , we should thus require 



Q (CO 



" (C^+Cmin) <0 



or 



Thus we have that 



Cky — Cky ^ Cn 



( hi + \q (cO 



(61) 



\Q\ 1/2 



dk x [5 (k x+ \Q\ 1 / 2 )+5{k x -\Q\ 1 / 2 )} 



V%Es a/— Cky — Cmin 

In the plane wave limit, we have then (converting to Ck y 

dCky 



(i) 

Q = 



8tt 2 C^ c 9/4 1 



, ,0^» \/Cl - Cn 



au(u 2 - ul) 2 



Ck y uJc9n z Cky + ujg' nz (Cky) ■ (62) 



We consider now two limiting cases, corresponding to weak static magnetic and electric 
fields respectively. 
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a. The weak magnetic field limit. The plane wave cross section in Eq. ()62j) can be 
simplified when the cyclotron frequency, u> c , is much smaller than the laser frequency, u, 
i.e., ijj c <C to. In this case, 



a 



(i) 



o 9/4 00 
3cr u uJc 



k 3 J2E 



E 



where we have defined 



a 



C J f 

a n,=0 ™m 



3ca; 3 



Ck v Cn 



(63) 



(64) 



in which cr° is the photodetachment cross section for H in the monochromatic field limit in 
the absence of any static fields, and k is the magnitude of the detached electron's momentum, 
k 2 = 2E f = 2{uj + Si). 

We note that our weak magnetic field result in Eq. (|63|) agrees with the formula of Peters 
and Delos (see Eqs. (3.6) and (3.7a) of Ref. |43j). Eq. (jo^j) agrees also with Fabrikant's 
result (see Eq. (53) of Ref. j^) except for the extra term in his formula that accounts for 
final-state interaction of the electron with the atomic residue. 

b. Weak static electric field limit. In the limit Es — > 0, we have that Cmin —oo. 
And we have also 



lim \f2E~ s J Ck y ~ Cmin = ^l /A \ 1 2uj r 



(£i + u) /lj c - (n z + -) 



(65) 



Substituting this result into Eq. (|62j) and carrying out the integration involving the Hermite 
polynomials, we obtain 

,2 



a 



(i) 

a=0,E s =0 



8tt 2 C?u; 2 c ^ r " 2 



CUj((jJ + U} c ) 2 



E 

n z =0 



(u - UJ C ) 



1 



2u c [(ei + u) /u c - (n z + |)] 



where the upper limit of summation, rii, is the largest integer that satisfies, 

Si + LU 1 

ni < 

[ uj c 2 

Eq. (|66p is exactly the same as Gao's result for the one-photon detachment cross section in 
a static uniform magnetic field (see Eq. (31) of Ref. 0|). 
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2. Cross Section for the Double Pulse Case 

From Eqs. (jUJ), P7f). (|5U|) and (J52j) . it is easy to show that for the double laser pulse 
case, the cross section is given by 



a 



(1 ) _ W 2 C?U} 2 C 
dbl „ ,(, ,2 , ,2^i 2 



cu;(u; 2 — u;;?) 2 

n z =0' 



; — ' poo 



X 



oo 



oo 



1 + COS [{S fi — Lj) T — f3]-= , 

^ l + cos^exp(-aVV2) ^ ^ - w) • (6?) 



We note that for r = /3 = 0, this formula reduces to the single pulse result in Eq. (j52J), as it 
should (cf. Eq. gTJ)). 



III. CONNECTIONS TO CLASSICAL CLOSED ORBITS 

In the previous section we have derived a general quantum mechanical expression for 
the (effective) photodetachment cross section for H~ by a short laser pulse in the presence 
of crossed static electric and magnetic fields. We have also shown that our plane wave 
limit result (given by Eq. (jH21)) reduces for the limiting cases of weak static magnetic (cf. 
Eq. (J63|) ) or weak static electric (cf. Eq. (J66|) ) fields to known results of others. Magnetic field 
strengths, B, that are readily available at present in the laboratory are weak in the sense that 
they satisfy the relation, u c uj. Therefore the quantum result for the photodetachment 
cross section in the plane wave limit given in Eq. (JBlIj) is of great interest owing to the 
possibility of experimental measurements with currently available technology. In this section 
we analyze this equation for the purpose of making connection with the classical closed orbits 
analyzed by Peters and Delos |42|. This connection will prove useful for interpreting some 
of the numerical predictions presented in the next section. 

For uj 3> uj c , the detached electron energy lies in the region of large n z . In this limit the 
integrand in Eq. (JBUj) becomes highly oscillatory, as may be seen by considering the large n z 
(Plancherel-Rotach) limit of the Hermite function, g nz 



where we have defined 



3: 

ok 1 : 



Uz + \ I (^^/l _ V 2 - arccos + ^ 



(68) 



r) = C fey /V2n~fl, (69) 
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and note that 

e ^ arccos t] ^ ir — e, e — > + , 

i.e., the argument of the Hermite function, g nz ((k y ), must lie between the classical turning 
points: 

-V2n z + 1 < Cfc, < v^wT+T- 



The function ^ (^Cfc y J that occurs in the integrand of Eq. (JUHj) may be calculated by differ- 
entiation of Eq. (jUHj) with respect to , as follows: 



\ (9?^ (9 / ~ \ 1 



x {| (1 - 7/T ! sin + (2n, + 1) (1 - r/ 2 )^ cos [5(n 2 ,r/)]} , (70) 

where we have defined the phase 

S(n z ,r)) = (™z + ~^ (yV l ~ V 2 ~ arccos r/) + (71) 

Assuming that n 2 3> 1, Eq. (fTU|) can be simplified (in particular, the first term within the 
curly brackets can be ignored in comparison with the second term), so that we obtain: 

9L (&,) * Jl(2n z )Hl-r) 2 )* C os[S(n z: r})]. (72) 

Owing to the fact that n z is large, the phase function S(n z ,r]) changes significantly as rj 
varies (cf. Eq. (fTTf ). so that g' nz [Ck^j oscillates rapidly as a function of (k y - From Eq. (}63|) we 
see that the magnitude of the photodetachment cross section will have the highest maxima 
when the squares of the various g' n _ (Xk v j functions that are summed (over n z ) have their 
maxima and minima in phase with each other, i.e., when neighboring phase functions differ 
by an integer multiple of 7r: 

S (n z ,r)(n z )) - S (n z - l,r](n z - 1)) ~ -f—S(n z ,r]) = jn, where j = 0, ±1,±2, ... 

dn z 

This condition is similar to that found for the largest local maxima in the photodetachment 
cross section of H~ in the presence of parallel static electric and magnetic fields j^. We 
compute the total derivative of S(n z ,r)) as: 

"■SM = ^ + a 4^- (73) 



dn z dn z drj dn z 
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The partial derivatives of the phase S(n z , rj) follow straightforwardly from the definition in 
Eq. (|71[). The partial derivative, drj/dn z , is calculated using the definition in Eq. (JBTJj) and 
the expression for Qk y obtained from the energy conservation condition, Efi — u = 0, together 
with Eqs. (I5fij) - (|58|) and Eq. (Jfilj) . After some straightforward algebra, one obtains: 

d 



dn. 



S{n z ,ri) = -arccosr/ + a/1 - 7] 2 ^->/2e^, 



J7T. 



(74) 



The condition (fTl^l for the highest local maxima in the photodetachment cross section (lfiT?|) 
may be re-written in terms of the scaled energy e, defined by 

2 



and the angle ip, defined by 



to obtain: 



COS (f + 



7T 



E s 



arccos rj, 



j + - I 7T = 0. 

2e V 2 ' 



(75) 



(76) 



(77) 



This result is identical to the classical equation expressing the relationship of the azimuthal 
angle (p and the scaled energy e for a closed orbit of an electron in crossed fields (see Eq. (3.12) 



of Ref. 



4*. 



The classical Hamiltonian corresponding to the quantum Hamiltonian ()12|) for the de- 
tached electron is given by 



H, 



els 



1 2 2 ( E,s\ 1 2 1 2 1 2 

-u c z + u c zlp y - — +-P z + - Px + -P y 



Denoting 



1 2 ^ E s 1-1 
2 P * + ~ Py + 2 Pz + 2 



£ *=2 P * + 2 



E s 

U C Z+ [Py 

UJ C 



2u>r 







UJ C Z + 


( p »-f). 







and introducing the following scaled coordinate, momentum, and time variables, 



(78) 



(79) 



(80) 



(81) 
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t = u c t, 



(82) 



Eq. ()78)1 may be rewritten as 



!~2 - 1 
^ = ^+2^+^" 2' 



(83) 



where E — ou^(u + e^/E'g is the scaled total energy and e is given by Eq. (|75|) in which the 
quantum energy e nz (cf. Eq. (J20*J) ) is replaced by the classical energy e z in Eq. (fTHJ) . 

With the help of Eqs. (|2U|) . (|19j) and ([bljh it is easy to show from the definition (J69j) that 



n 



Py ~ E s/ 



0J r 



(84) 



which can be rewritten, in terms of scaled energies, as 

e-{E- 1/2) 

v = — 



(85) 

by using the energy conservation equation (|5S|l and the fact that p x = for closed orbits. 
Substituting Eq. into Eq. (fTUj) . we can rewrite Eq. (|77j) as 

■e-(S- 1/2)" 



A(e) = y 2e — [e — (E — 1/2)} Z - arccos 



2e 



J7T. 



(86) 



For a given scaled total energy E, the number of the solutions of Eq. (J86|) gives the total 
number of closed orbits. The return time of a closed orbit in crossed fields is given by 

T ret = 2(^ c )- 1 v / 27cos< A (87) 

which can be rewritten with the help of Eqs. (f?T)|) and (f83J) as 



T ret = 2(lu c )~ 1 J2e -[s-(E- 1/2)]' 



(88) 



As discussed in |42j, there exists a very important group of closed orbits whose total 
energies are given approximately (in the large energy limit) by 



TT^Ej 

2 ul 



E)^ 



j + 



71- 



(89) 



where j = 1,2,3, .... These are called boundary energies, because for each j a new closed 
orbit appears at the energy given by Eq. ()89|) and for higher total energies this newborn 
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closed orbit will split (or "bifurcate" ) into a pair of closed orbits with two different energies 
and return times, given by Eqs. (j8T?j) and ([88)1 . 

Actually, each boundary energy defines the onset of large oscillations in the cross section. 
However, the largest amplitude oscillation in the cross section occurs at a slightly higher 
energy at which a different type of closed orbit occurs that has a truly circular motion in 
the drift frame in the y-z plane. The energy of this orbit may be obtained by setting the 
initial momentum along the y axis equal to the drift velocity, i.e., 

Pi = -• 00) 



From the energy conservation equation ()83|) and the fact that p x = for a closed orbit, we 
have e = E — 1/2. Substituting this result into Eq. (JHUj) gives 

V2~e=(j + ^)n, (91) 

which in unsealed variables corresponds to a total energy equal to 



— TT 2 E 2 g 



3 + 



1 

+ — 



(92) 



Comparing Eqs. (|8*9*j) and (J9*2|) . one sees that the energy difference between the boundary 
orbits and the orbits having p® equal to the drift velocity is 2Eg/u 2 , independent of the value 
of j. Boundary closed orbits satisfy dA(e)/de = 0, which gives the relationship e = E + | 
in the large energy limit. Closed orbits for which Eq. (|90|) applies have s = E — |. From 
Eqs. (j88|) and (J9T|) . we find for these latter orbits that 

T ret = (3 + l)T B , (93) 

where Tb = 2h/oj c is the cyclotron period and j is a positive integer. This formula is very 
similar to that obtained for the case of parallel static electric and magnetic fields [3, 0| , 
in which the largest oscillation amplitude of the cross section corresponds to classical orbits 
for which for an electron is ejected along the static field direction and reflected by the static 
electric field such that its return time satisfies T ret = jTg. In the parallel fields case, the 
motion in the plane perpendicular to the magnetic field is simply cyclotron motion with 
period Tb- Classical closed orbits having a return time equal to an integer multiple of Tb 
are associated with the largest oscillations in the cross section 
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FIG. 2: Photodetachment cross section for H~ in the plane wave limit for a static magnetic field 
B = 1 T and three different values of the static electric field strength. Results are plotted versus 
detached electron kinetic energy above the zero field detachment threshold, uJ + Ei, up to 2.7 cm -1 . 



In the crossed fields case, however, the situation is much more complicated. However, 
since e z (cf. Eq. ()79|0 is conserved, when the detached electron has an initial momentum 
Py given by Eq. (JiJUj) (and an initial position z — 0), the initial momentum along the z 
axis takes its maximum value. This implies that this particular closed orbit starts out (in 
the drift frame) aligned with the laser polarization direction. These circular orbits (in the 
drift frame), having energies given by Eq. (}9"2*j) . are associated with the largest amplitude 
oscillation of the cross section. 



IV. RESULTS AND DISCUSSION 

In this section we present numerical results based on the quantum mechanical theoretical 
formulation presented above. We present first plane wave limit results for the photodetach- 
ment cross section of H~ in the presence of crossed static electric and magnetic fields over 
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FIG. 3: Same as Fig. [21 but for three higher static electric field strengths and detached electron 
kinetic energies up to 5 cm -1 . 



a much larger energy range than in prior works 



Q, 42, Q 



This large range allows us to 



421, both in 



demonstrate very clearly the signatures of the predicted classical closed orbits 
the energy spectrum and in the time (i.e., Fourier transform) spectrum. We examine next 
the short laser pulse case, demonstrating first the effects of laser pulse duration on the pho- 
todetachment cross section. We then examine the detached electron wave packet dynamics 
in the y-z plane and the possibility of modulating the detachment cross section by pump 
probe (Ramsey interference) techniques. The connection between the time development of 
the quantum wave packet of the detached electron and the predicted classical closed orbits 
is also discussed. 
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A. Photodetachment Cross Section in the Plane Wave Limit 

1. Static Electric Field Dependence for Near Threshold Energies 

In Figs. El and 01 we present the photodetachment cross section for H~ for a static magnetic 
field, B — 1 T, and six different values of the static electric field, Eg- Our quantum theory 
predictions are obtained from the plane wave limit result given in Eq. ()62|). One sees in Fig. EI 



that, as noted by Fabrikant 4l(, even a very small static electric field removes the known 
singularity in the detachment cross section for energies corresponding to integer multiples 
of the cyclotron frequency in the pure magnetic field case (see, e.g., |5j|). In particular, for 
Eg = 0.5 y/cm, the behavior of the cross section is very similar to that of the pure magnetic 
field case |57( (to which our results reduce in the limit of zero static electric field, as shown in 
Sec. II. F. 1(b) above), but without the cyclotron singularities. On the other hand, beginning 
with Eg = 7 V/cm, the oscillatory modulation of the cross section by the static electric field 
becomes obvious. As shown in Fig.EJthe frequency of this modulation decreases as the static 
electric field magnitude increases, just as is found for the case of a pure static electric field 
or for the case of parallel static magnetic and electric fields (see, e.g., 0]). One sees also 
that the cross section becomes non-zero at the zero-field threshold owing to the lowering of 
the threshold by the static electric field. For the present crossed static magnetic and electric 
field case, the modulation of the cross section becomes increasingly complex the higher the 
maximum total energy Ef becomes. For a maximum detached electron kinetic energy of 30 
cm -1 , Fig. |U shows that the oscillatory modulations differ above and below approximately 
15 cm -1 . For energies below 15 cm' 1 , there exists only a sinusoidal modulation. Above 
about 15 cm" 1 , the modulation consists of more than one frequency and becomes more 
complicated the higher in energy one looks. 

In order to examine these structures in detail, it is instructive to plot only the oscillatory 
part of the cross section, which is defined by 

a osc = aE - a°, (94) 

where cr£^ is the total cross section in the plane wave limit (given by Eq. (jfi2j) ) and a is the 
photodetachment cross section in the absence of any external static fields (given by Eq. (jHljl ). 
Figure El shows the oscillatory part of the cross section over three different energy ranges, 
corresponding to total energies up to 60 cm -1 , 180 cm -1 , and 500 cm -1 respectively. While 
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FIG. 4: Photodetachment cross section for H in the plane wave limit for B = 1 T, .Eg = 60 
V/cm, and detached electron kinetic energies up to 30 cm -1 . 

the oscillatory modulations of the cross section become increasingly dense and complex 
as the total energy increases, we see also that clear patterns in the spectra emerge and 
become more obvious the higher in energy we look. The onset of these repetitive patterns 
is indicated in each panel by the vertical dotted lines, which represent the locations of the 
boundary energies [42j defined by Eq. The peak amplitudes are indicated by the open 
triangles at the energies defined by Eq. which correspond to the locations of circular 
classical orbits in the drift frame. The connection of these classical closed orbits and our 
quantum mechanical cross sections can be most easily investigated in the time domain, which 
we consider next. 



2. Fourier Transform Spectra and Closed Classical Orbits 

The Fourier transform of the oscillatory part of the photodetachment cross section, a osc , 
is presented in Fig. |H1 for increasing values of the maximum total detached electron kinetic 
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FIG. 5: The oscillatory part of the cross section, cr osc (cf. Eq. (EH)), for B = 1 T and E s = 60 
V/cm for u! + e% ranging from (a) to 60 cm , (b) to 180 cm -1 ; (c) to 500 cm -1 . Dashed 
lines indicate the boundary energies (cf. Eq. (|89|)) at which a new closed orbit appears [3] and the 
open triangles indicate the energies at which the amplitude of the oscillatory part of the quantum 
cross section is expected to have a local maximum (cf. Eq. (|92[)). 

energy, _E™ a:r . Times are given in units of the cyclotron period, Tg. We see from Fig. E(a) 
that at the lowest maximum energy, only one peak appears in the time spectrum. The 
peak position indicates the return time (i.e., orbit period) of a closed classical orbit having 
an energy of 10.98 cm -1 . As E™ ax increases to 15.91 cm -1 in (b), another peak emerges. 
This corresponds to the first classical boundary energy (cf. Eq. (|89Jl ) for j — 1. Unlike 
the case of classical dynamics, however, where the boundary energy is sharply defined, our 
calculations indicate that the second peak in Fig. E(b) begins to appear around the energy 
13.5 cm -1 . Note also that the first peak in (b) shifts to the right as compared to that in 
(a). This indicates that the return time of the first closed orbit increases when the total 
energy increases. In panel (c) we observe that the second peak increases in magnitude while 
the first peak decreases in magnitude. In panel (d), when the total energy equals 21.97 
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FIG. 6: Fourier transform spectra for the oscillatory part of the cross section (cf. Eq. (j94[) ) for 
different maximum energies, E™' 1 * = u> + £j, as indicated in each panel. See text for a detailed 
description. 

cm -1 , one observes the bifurcation of the second peak that first appeared in panel (b). For 
j^max _ 43 4g cm _1 ; in panel (e), one sees the appearance of a third peak around 2.4 Tg 
and notice that the width of the splitting of the second peak increases as compared to that 
in (d). The left and right boundaries of the split and broadened second peak correspond to 
the return times of the two bifurcated orbits at the maximum available total energy EJ 1 ^. 
The serrated U-shaped region between the left and right boundaries of the split second peak 
correspond to the return times of the two bifurcated orbits for lower total energies (e.g., 
such as those two shown in panel (d) for a total energy of 21.97 cm -1 ). Finally, in panel (f) 
for a slightly higher total energy we observe that the third peak grows in magnitude relative 
to the first and second (split) peaks. 

Consider now a much larger total final state energy, EJ 1 ^ = 500 cm -1 . The oscillatory 
part of the cross section is shown in Fig. Efc). In this figure, the dashed lines correspond 
to the seven boundary orbit energies, given by Eq. (jBHJ), that appear for energies up to 500 
cm -1 , and the triangles correspond to the seven circular orbit (in the drift frame) energies, 
given by Eq. (}92|) . The Fourier transform spectrum for energies in the range from - 500 
cm -1 is shown in Fig. Ufa), while the Fourier transform of only the part of the spectrum in 
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TABLE I: Numerical solutions of Eq. (|86|) for the energies, e- (in units of cm" 1 ), for the closed 
orbits that exist for a total energy u>+ei = 500 cm -1 (302.869 in scaled units) and the corresponding 
closed orbit periods, Tp (in units of Tg), calculated from Eq. (|88[). Note that there is only one 
solution for the case j = 0. 
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461.459 


462.389 


463.987 


466.337 


469.595 
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480.398 


491.231 


T i 


0.959 


1.918 


2.876 


3.831 


4.783 


5.730 


6.665 


7.571 
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542.137 


541.0371 


539.127 


536.264 


532.160 


526.143 


515.609 






1.041 


2.082 


3.126 


4.172 


5.224 


6.286 


7.378 



the range from 400 - 500 cm -1 is shown in Fig. [7(b). In (a), the open circles denote the 
return times (periods) of the 15 closed classical orbit solutions of Eq. ()86|) that exist for a 
total energy of 500 cm" 1 . These periods are calculated using Eq. (}88|) and the results are 
given in Table |J together with the corresponding orbit energies. (Note that for each j > 
and for a total energy E not equal to one of the boundary energies given by Eq. (jEHJ), Eq. (jHHJ) 
has two solutions.) The open triangles, on the other hand, correspond to the circular orbits 
in the drift frame corresponding to the local maximum amplitudes of the oscillatory part of 
the cross section; their return times are given by the very simple Eq. (j§3l . 

Note that the bowl-like structures appearing in Fig. UJa) above each open triangle result 
from the fact that closed orbits have different return times for each different total energy and 
from the fact that the Fourier transform spectrum in this figure results from a large range 
of total energies, i.e., from to 500 cm -1 . When we calculate the Fourier transform of the 
oscillatory part of the cross section over only the limited energy range, 400 cm -1 < uo + Ei < 
500 cm -1 , as in Fig. UJb), then we observe that the first 13 peaks are approximately located 
at the positions of the first 13 closed classical orbit periods given in Table HI which were 
calculated for a maximum total energy of 500 cm" 1 . The energy region < uo + £, < 400 
cm -1 is thus inferred to be responsible for the bowl-like structures in Fig. Ufa) owing to 
the shifts of the lowest 13 classical orbit periods (having j < 6) for lower total energies. 
The bowl-like structure remains between the 14th and 15th closed classical orbits (having 
j = 7) as this pair of orbits first occurs above a total energy of approximately 455 cm" 1 (cf. 
Fig. Efc)). We also observe from the data in Figs. Eland UJ that the oscillation amplitude of 
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FIG. 7: Fourier transform spectra of the oscillatory part of the cross section, <7 OS c (cf. Eq. (|94[)). 
given in Fig. 03 (c) calculated over two different energy ranges: (a) < uj + £i < 500 cm" 1 ; (b) 400 
cm -1 < uj + £i < 500 cm" 1 . In (a) the open circles indicate the return times (i.e., orbit periods) of 
the 15 closed orbits for uj + Si = 500 cm -1 (see TableQ]). Also in (a), the open triangles indicate the 
return times of the circular closed orbits (in the drift frame) having p® = I (cf. Eqs. (|81|) and (|90ll ): 
these return times are given by Eq. (|93|) . 

the cross section becomes larger as the total energy increases. 



B. Detachment by Short Laser Pulses 

Photodetachment by means of one or more short laser pulses differs from that by a 
monochromatic laser. Most obviously, the pulse bandwidth affects the measured spectrum 
of detached electrons. In addition, short laser pulses produce localized detached electron 
wave packets whose motion in crossed fields can be investigated and compared to classical 
predictions 0, Q, Q. Most interesting, perhaps, is the possibility of controlling the 
modulation of the detachment spectrum by variation of the parameters of one or more laser 
pulses. In the rest of this section, we examine each of these topics in turn for the crossed 
static electric and magnetic field case. 

We note first, however, several previous works on related problems. Ramsey interference 
effects resulting from photodetachment of H~ by two short, coherent laser pulses as a function 
of their relative phase was examined by Wang and Starace for the case of a static electric 
Md Q aad the case of paraM static electric aad magaetic he.ds Q. The latter w 0rk 
showed that a large modulation of the effective detachment probability can be achieved by 
optimizing the static field magnitudes and the time delay between laser pulses, as follows: 
the field magnitudes should be such that the classical time for reflection of an electron back 
to the origin by the static electric field equals an integer multiple of the harmonic oscillator 
period for electron motion in the static magnetic field; also, the time delay of the second 
pulse should coincide with the classical time for the electron's return to the origin. For 

n 

the case of a single short laser pulse, Du 62] examined the photodetachment of H in the 
presence of a static electric field using modified closed orbit formulas. He showed that when 
the laser pulse duration is shorter than particular closed orbit periods, then those orbits no 
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FIG. 8: Effective photodetachment cross section of H~ (cf. Eq. (|53|l) by a single laser pulse of the 
form (jBJ) with four different pulse durations (cf. Eq. (J7J) in the presence of crossed static electric 
and magnetic fields, Eg = 15 V/cm and B = 1 T. Results are plotted vs. electron kinetic energy 
beginning from the zero-field threshold. Also shown (dotted line) is the photodetachment cross 
section for the case of a continuous (monochromatic) laser without any external static fields present. 



longer contribute to the photodetachment spectrum. Finally, Zhao et al. [63j have derived 
a uniform semiclassical formula for the photodetachment cross section of a negative ion by 
a short laser pulse for the case of parallel static electric and magnetic fields. 



1. Pulse Duration Effects 

The fundamental difference between using a short laser pulse and using a continuous 
(monochromatic) laser is the bandwidth of the short laser pulse. In the former case, the 
laser pulse will excite a group of final states that form an electron wave packet, whereas 
in the latter case only a well-defined final state will be reached. In our present case of 
detachment in the presence of crossed static electric and magnetic fields, the spacing of 
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Landau levels is very small (0.93 cm" 1 for B = 1 T), so that one expects that even a quite 
long pulse having a duration of several picoseconds will have considerable finite bandwidth 
effects on the detached electron wave packet and its dynamics. 

In Fig. El we present the effective total cross section (cf. Eq. (|53jl ) for a laser pulse of the 
form (jHJl (with r = j3 = 0) and four different pulse durations (cf. Eq. (JZj)) in the presence 
of crossed static electric and magnetic fields of magnitudes Eg = 15 V/cm and B = 1 T. 
As we see from Fig. |HJ for the longest pulse duration, 240 ps, the effective cross section is 
identical with that for a (monochromatic) plane wave. As the pulse duration decreases, the 
modulation of the cross section is suppressed, beginning at the highest energies shown and 
progressing to structures at lower energies. Thus, when the pulse duration is reduced to 45 
ps, the modulation structure beyond the energy 2.3 cm -1 is largely suppressed. As the pulse 
duration is further reduced to 30 ps, even the modulation between 1 and 2 cm -1 decreases 
in magnitude. At the shortest pulse duration, 15 ps, the oscillatory structure completely 
disappears and the effective cross section becomes a smooth curve passing through the 
oscillatory cross sections for longer pulse durations. In this case, the cross section is nearly 
identical to the one for detachment by a continuous (monochromatic) laser in the absence 
of any external static fields, as shown by the dotted line in Fig. |HJ The major difference 
between these cases occurs near the threshold: our short pulse effective cross section is finite 
at the zero static field threshold, while the monochromatic field cross section, in accordance 
with Wigner's threshold law, is zero at threshold. This difference is due to the lowering of 
the detachment threshold by the static electric field. 

It is interesting to relate these changes in the structure of the effective detachment cross 
sections as a function of pulse duration to the energy positions of the known classical closed 



orbits 



43J. For the maximum energy 2.7 cm 1 considered here, there are three closed 



orbits available for the static field parameters we employ. These three orbits have return 
times (periods) of 30.69, 40.94 and 60.72 ps. As the energy decreases to 1.8 cm -1 , the 
return times decrease to 28.94, 43.32 and 56.19 ps respectively. As the energy is further 
decreased to 1 cm -1 , there are only two closed orbits having return times of 26.95 and 48.58 
ps respectively. We observe in Fig. [H] that as the laser pulse durations become shorter than 
the closed orbit periods, the structure of the effective cross section is reduced. In particular, 
for the shortest pulse duration, 15 ps, which is smaller than any closed orbit period, all 
structure has disappeared. 
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FIG. 9: (color online). Contour plot snapshots of detached electron wave packet motion in the y-z 
plane for the case of zero initial momentum in the x direction [cf. Eq. (|95|)]. The static electric 
and magnetic field strengths are 60 V/cm and 1 T respectively. The laser pulse duration is T p = 2 
ps and t = corresponds to the end of this pulse. The total electron energy is Ef = 8 cm -1 . 

2. Wave packet dynamics 

We examine here the dynamics of a detached electron wave packet produced by a short 
laser pulse under the influence of crossed static electric and magnetic fields. As discussed 



in Refs. |42i |43j, the oscillatory part of the photodetachment cross section produced by a 
monochromatic laser field may be associated with those closed orbits that exist for a given 
value of the photon energy. As discussed in the previous section, the oscillatory part of the 
effective cross section produced by a short laser pulse is suppressed when the pulse duration 
is smaller than the classical orbit periods of those orbits that exist at the energy being 
considered. (See for the related pure static electric field case.) From the correspondence 
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FIG. 10: (color online). Same as Fig.|§]but for a total detached electron energy of Ef = 15.9 cm . 

between classical and quantum mechanics, we expect to observe that the detached electron 
wave packets produced by short laser pulses will trace the paths of allowed classical closed 
orbits. 

In order to illustrate the quantum wave packet motion corresponding to the classical dy- 
namics, we shall only consider two-dimensional quantum motion by imposing the restriction 
p x = 0. As discussed above, the x component of the momentum has to be zero in order 
for there to be any closed orbits. In this case, the time-dependent electron wave packet in 
coordinate space is given by 

i r°° 

^(0,y,z,t) = - ^ k (0,p y ,z,t)e^y, (95) 

which is the Fourier transform of Eq. (}4T)|) . taking p x = 0. A similar Fourier transform can 
be employed for the double pulse case in Eq. (jlSJ). 

In Figs. I9TITT1 we present snapshots of the detached electron wave packet for the case of 
a single laser pulse of the form © with pulse duration T p = 2 ps. The time evolution 
starts from to = —T p . The static electric and magnetic field strengths are taken to be 60 
V/cm and 1 T respectively in all cases. Note that the cyclotron period is Tb = 35.72 ps for 
B = 1 T. In Fig. El we take the total energy Ef to be 8 cm -1 . There is only one classical 
closed orbit for this energy and these static field parameters. The return time of the closed 
orbit is calculated to be T® et = 24.1 ps (0.674 Tb). In Fig. HJa), we see that two electron 
wave packets are created at the peak intensity of the laser pulse on either side of the z = 
axis, which correspond to electrons being ejected either along or opposite to the direction 
of increasing static electric field, E$ (cf. Fig. HJ. After the end of the pulse in (b), the two 
electron wave packets move apart. However, as time increases we see in (c) that both wave 
packets are turned back by the external magnetic field. As shown in (c) and (d), each wave 
packet undergoes considerable spatial spreading. The most interesting plot is shown in (e), 
where we see a large portion of the left hand wave packet sweep through the residual core 
(at the origin). We note that the time corresponding to this snapshot is exactly the return 
time of the only classical closed orbit in the present case. It is return of this piece of the 
wave packet that leads to the regular sinusoidal oscillation one sees in Fig. Efa) below 15.9 
cm -1 . As the time approaches 1 T B in (f), we see that the two wave packets refocus on the 
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FIG. 11: (color online). Same as Figs. lTuT e^. (f), (g) and (h) except that the wave packet amplitudes 
are shown here in three-dimensions rather than as contour plots. Note the scale change in panel 
(b), which shows the re-focusing of the wave packet amplitude along the drift axis (i.e., away from 
the atomic core at the origin). 

positive y axis (the direction of drift motion) and that they pass through each other and 
continue their rotational motion during the next cyclotron period. However, owing to the 
drift motion along the y-axis, we see in (h) that the left hand wave packet is no longer able 
to return to the atomic core in the second cyclotron cycle for this total energy. The two 
wave packets do refocus further along the positive y-axis again at 2 Tg, as shown in (i). 

We note that the refocusing and the drift of the electron wave packets are exactly anal- 
ogous to the classical dynamics discussed by Peters and Delos in (3]. They showed that 
classical orbits with different initial conditions will refocus at various points along the drift 
axis. We note also that pump-probe experimental studies of the related problem of the 
motion of pump-laser-produced Rydberg-state electron wave packets in the Rubidium atom 
in the presence of crossed fields have found enhancements of probe laser-produced ionization 
signals when the delay of the probe laser equals the orbital period of the appropriate closed 
classical orbit for this related problem Q]. 

We present similar snapshots in time for the increased energy of 15.9 cm -1 in Fig. E3 At 
this energy, the first so-called boundary orbit j^J may be populated (cf. Eq. (jHHJ) for j = 1). 
There are thus two classical closed orbits that exist at this total energy, whose return times 
are 27.4 ps (0.766 Tg) and 48.8 ps (1.366 Tg). In Figs. HUf e) and lTOT g) we see that different 
parts of the quantum electron wave packet return to the atomic core at these two times. The 
most distinctive feature of the part of the electron wave packet that returns to the origin 
at approximately 49 ps (cf. Fig. ITOT g)) (which corresponds to the higher energy, classical 
boundary orbit) is that most of the arc in (g) passes through the atomic core at the origin 
(which we have confirmed by observing the motion of the electron wave packet on a finer 
time scale). We note also that the energy of the classical boundary orbit corresponds to the 
abruptly increased amplitude of the oscillatory part of the cross section seen in Fig. Efa) 
around the energy location 15.9 cm -1 indicated by the first dashed line. 

The fact that electron wave packet amplitudes return to the region of the atomic core 
implies the possibility of modulating the detachment cross section, analogously to the case 



35 



of a monochromatic laser, as shown in Figs. H] and |3J However, the present wave packet 
studies show why for the crossed field case the modulation of the cross section is very small. 
Consider the three dimensional wave packet snapshots in Fig. ^2 which are calculated for 
times corresponding to those in Figs. ITUT e). (f), (g), and (h). Owing to the spreading of the 
electron wave packet, when it returns to the origin the part of the wave packet that overlaps 
the origin is nearly two orders of magnitude smaller than the probability in (f), in which 
the wave packet re- focuses along the drift axis, i.e., away from the origin. Because of such 
wave packet spreading and drift away from the origin, modulation of the photodetachment 
cross section in the crossed static field case is necessarily small. For similar reasons, the use 
of short pulse, pump-probe type techniques to control the photodetachment cross section in 
the crossed static field case also results in only small modulations of the cross section, as we 
discuss next. 



3. Pump-Probe Coherent Control of the Effective Photodetachment Cross Section Using Short 

Laser Pulses 

The idea of using laser pulses shorter than electron wave packet orbital periods to control 
electron wave packet motion was initially formulated theoretically for Rydberg (i.e., bound) 
electron wave packets This idea was extended theoretically to photodetached (i.e., 



external static fields, including static 
fields. Experimentally, short pulse, 



continuum) electron wave packets in the presence o 
electric |6jJ and parallel static electric and magnetic 
pump-probe studies of photodetachment of 0~ in the presence of a static magnetic field 
demonstrated Ramsey interference between photodetached electron wave packets j^. Such 
Ramsey interference may also be demonstrated in the present crossed electric and magnetic 
field case. 

In Fig.^Jwe present the effective total photodetachment cross section (cf. Eq. (|67)l) as a 
function of the relative phase (3 and the time delay r between two laser pulses (cf. Eq (J41|) ) 
for a total detached electron energy of 15.9 cm -1 and for E$ = 60 V/cm and B — 1 T. 
The pulse duration T p of both pulses is taken to be 4 ps. Fig. H2Ta) shows the dependence 
of the effective cross section on the relative phase, /3, for six time delays, r, between the 
pulses. One observes that the modulations of the cross section have local maxima for time 
delays of 0.77 T B (~ 27.5 ps) and 1.37 T B (~ 49 ps), which are precisely the return times 
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FIG. 12: The effective cross section for the double laser pulse case as modulated by: (a) the relative 
phase between the two pulses for several time delays, as indicated; (b) the time delay between the 
two pulses for two fixed relative phases, and ix. 

of the two allowed classical closed orbits for a total electron energy of 15.9 cm -1 . However, 
the modulation of the cross section for the larger time delay is much greater than for the 
smaller time delay, which is consistent with the extent of electron wave packet overlap with 
the origin shown in Figs. IToT e) and (g). In other words, in the latter case a large portion 
of the wave packet passes over the origin, which makes the Ramsey interference with the 
newly produced wave packet amplitude (due to the second pulse) of greater amplitude. 

Fig- E(b) shows the dependence of the effective cross section on the time delay, r, for 
two relative phases, (3, between the pulses: and 1 n. Fig. [T2T b) clearly shows that the 
maxima and minima in the effective cross section as a function of the time delay between 
the pulses occur for time delays of 0.77 Tb (~ 27.5 ps) and 1.37 Tb (~ 49 ps), which are 
the orbital periods of the two allowed classical closed orbits. We see once again that the 
modulation of the effective cross section is much larger for the classical closed orbit having 
the larger time delay, as explained above. 
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V. CONCLUSIONS 

In this paper we have presented a detailed quantum mechanical analysis of detachment 
of a weakly bound electron by a short laser pulse in the presence of crossed static electric 
and magnetic fields. For specificity, we have chosen the parameters of the initial state of 
the weakly bound electron as those appropriate for the outer electron in H~. In particular 
we have presented an analytic expression for the final state electron wave function, i.e., the 
wave function for an electron moving in the field of a laser pulse of arbitrary intensity as 
well as in crossed static electric and magnetic fields of arbitrary strengths. The general 
detachment probability formulas we present may therefore be used to analyze multiphoton 
detachment in crossed fields (although we have not presented this analysis here, but instead 
have focused on the weak laser field case). 

Based upon our analytic results for the detachment probability by a short laser pulse, 
we have defined an effective detachment cross section for the short pulse case that is shown 
to reduce, in the long pulse limit to results of others for the monochromatic, plane wave 
case. Our effective cross section formula allows us to demonstrate the effects of the laser 
pulse duration, such as, e.g., that for pulse durations shorter than the period of a particular 
classical closed orbit, the features of that closed orbit in the photodetachment spectrum (for 
the plane wave case) will simply vanish. By means of a stationary phase analysis, we have 
derived a condition for the existence of closed classical orbits that agrees exactly with that 
obtained by Peters and Delos by a purely classical analysis j^j]. We have also illustrated the 
bifurcation of the closed classical orbits at the so-called boundary energies j^j] by Fourier 
transforming the oscillatory part of our quantum cross section (in the long pulse limit) over 
various ranges of the final state electron energy. 

Finally, our analysis of the motion of detached electron wave packets produced by a 
short laser pulse provides a direct comparison of quantum and classical features for the 
crossed static electric and magnetic field problem. We find that the dynamics of our two- 
dimensional detached electron wave packets are consistent with the predictions of closed 
classical orbit theory [3] . We have also shown that wave packet spreading and the fact that 
wave packet refocusing only occurs at the origin in the drift frame means that control of 
electron detachment in crossed static fields by means of laser pulses is less effective than in 
the parallel static electric and magnetic field case 
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APPENDIX A: ANALYTICAL WAVE FUNCTION FOR FREE ELECTRON 
MOTION IN A LASER FIELD AND CROSSED STATIC ELECTRIC AND 

MAGNETIC FIELDS 



In this appendix we give the details of the solution of Eq. (|10[). which describes free 
electron motion in a laser field in the presence of crossed static electric and magnetic fields. 
The configuration of the external fields is shown in Fig. ^ The solution is clearly separable 
in momentum space and has the following form: 



"0/ (P,*) = tPx(Px,t)lpy(Py,t)lp z (p z ,t)eXp 

= 5{p x - k x )5(p y - k y )ip z {p zi t) 



x exp 



(Al) 



where the z component of the final state wave function satisfies the following equation: 



1 2 d 2 . ( E s \ d 1 2 1 . . 



■ d , , , 
i-Qj.iPzKPzit) 



In order to solve Eq. (1A2I) . we make the following substitution: 



ipz{p z ,t) = exp 



-i—(k y - — )p« 

UJc UJ r . 



MPz,t). (A2) 



(A3) 



which serves to eliminate the term involving the first derivative of p z from the equation 
satisfied by i> z i{p z ,t), 



■ d , , , 

^Wzl{Pz,t) 



2" c dp 2 z ' 2 
Upon making the substitution, 

i/>zi(Pz,t) = exp 



1 2 d 2 1 o 1 , , , 1 , 

-UJ 2 ^ + -Pi + -p z A L (t) - - ( ky 



UJ C 



4>zl(Pz,t)- 



(A4) 
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2 
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(A5) 



Eq. (|A4J) gives the following equation for ip Z 2{pz,t)'- 



■ d , , , 

l^z2{Pzit) 



1 2 92 1 2 1 „ / N 



^z2{Pz,t). 



(A6) 



Eq. (IA6I) has the form of the equation for a forced harmonic oscillator, which can be 

n 

solved exactly [67J: 



4>z2{p z ,t) = exp 



L(t')dt' 



-1/4 



(V2( Pz 



(A7) 
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where ( Pz and e nz are given by Eqs. (fTHj) and (plj) respectively. The function is related to 
the vector potential of the laser pulse and satisfies the differential equation given in Eq. ()22|) . 
The function L(t) is defined in Eq. (|21j) in terms of £,(t). In Eq. (|A7|) . the function g Hz (x) 
is defined by Eq. (JHJ). Note that g n X x ) * s normalized, i.e., 



/oo 
[g nz fdx = 1. 
-oo 



(A8) 



We note also that Qn z (x) is simply proportional to a parabolic cylinder function (see 
Eq. 19.13.1 of Ref. |^|): 



9n z {x) = i , U \-n z - -V2x 
Combining Eqs. (|A3j) . (JA5|) . and (|A7jl . we find for the solution of Eq. fjA"2|) . 



(A9) 



i>z(Pz,t) = exp -ib(k y ,t)V2( Pz u c 1,A g nz (y2C, Pz 



x exp 



- i^Cfe^(0 + i\u£lt + i / L(t')dt' 



(AlO) 



where b(k y ,t) and Cfc H are given by Eq. (|T7|) and (fTU|) respectively. Finally, substituting 
Eq. (jA10|) into Eq. ()A1|) gives us the analytical expression for the final state wave function 
in momentum space given in Eq. f)13j) . 



APPENDIX B: SOLUTION OF EQ. FOR f (t) 



In this appendix we present the solution for the laser- field-dependent function £(t), which 
satisfies the differential equation (j22J). The corresponding Green's function satisfies, 



G(t,t') +uj 2 c G(t,t') = 5{t-t'), 



whose solution is: 



G(t,t') 



2tt 



oo - 



du'. 



In terms of this Green's function, the solution of Eq. (J22|) is: 
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-Mt') 

c 
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(B3) 
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Using the definitions in Eqs. © and Q, the integral over t' can be evaluated (using inte- 
gration by parts) to obtain, 



dt' 
n E 



~A L (t>) 
c 



-iui't' 



icu' 2ia 

in which the function W{uj') is defined by 



and where it has been assumed that A L (— oo) = A L (oo) = 0. 
Thus £(t) is now given by 
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where we have defined the following three integrals over a/: 
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1. Calculation of I (t), h(t) and I 2 (t) 



Consider first h(t), which satisfies the following relation: 
d 



dV 



du'w^'y^'-^ 1 



i 
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Thus, we have that 
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Replacing uj c by — u c in the above formula, we get 



h{t) 
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Setting u; c equal to zero in Eq. (|B10J) and multiplying by a factor of —2, we get 



+ 



erf ( a(t _ T) + g + 1 



(B12) 



2. Some Relations for and Its Derivatives 



From Eqs. (jB5|) and ()B9|) (and equations similar to (jB9|) for io(^) an d hit)), we have that 



4ai 

E UJ C ^ 

4a 



I (t) + h{t) + I 2 (t) 
[h(t) - h(t)} . 



(B13) 
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The second derivative of is given by 



4a 
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Writing now 



where 



i=0 



Eqs. flBUjl and jEHj) give 



e(0 = ^ c [^(t)_e 2 (t)] 
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Substituting these equations into Eq. (j22j) gives: 

= — A L (t). 
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3. Long Pulse Approximation for 



In the long pulse case (in which a /to <C 1), s imp lified expressions can be obtained for the 
integrals Ii(t) in a way similar that used in Ref. For example, starting from Eq. (|B9)1 one 
can expand the sine function in terms of exponentials and then do the resulting integration 
by parts, dropping terms that are of order a/oo or higher, as follows: 

ft 
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In a similar way one may obtain the following approximate expressions for ^(t) and Io{t)'- 



h{t) 
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LO + UJ C 
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LO — L0 r 



Aia 1 _ a 2 (t _ T) 2 . 
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(B23) 
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Substituting these long pulse approximations for the integrals into Eq. (jB5|l . the 
function in the long pulse limit is then given by Eq. 
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